Unusual features of coarsening with detachment rates decreasing with cluster mass 
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(Dated: December 17, 2008) 

We study conserved one-dimensional models of particle diffusion, attachment and detachment 
from clusters, where the detachment rates decrease with increasing cluster size as 7(171) ~ m~ k , 
k > 0. Heuristic scaling arguments based on random walk properties show that the typical cluster 
size scales as (t/lnt) z , with z = l/(fe + 2). The coarsening of neighboring clusters is characterized 
by initial symmetric flux of particles between them followed by an effectively assymmetric flux due to 
the unbalanced detachement rates, which leads to the above logarithmic corrections. Small clusters 
have densities of order ( _mz ( 1 ) j with z (1) = k/ (k + 2). Thus, for k < 1, the small clusters (mass 
of order unity) are statistically dominant and the average cluster size does not scale as the size of 
typically large clusters does. We also solve the Master equation of the model under an independent 
interval approximation, which yields cluster distributions and exponent relations and gives the 
correct dominant coarsening exponent after suitable changes to incorporate effects of correlations. 
The coarsening of typical large clusters is described by the distribution P t (m) ~ l/t y f (m/t z ), with 
y = 2z. All results are confirmed by simulation, which also illustrates the unusual features of cluster 
size distributions, with a power law decay for small masses and a negatively skewed peak in the 
scaling region. The detachment rates considered here can apply in the presence of strong attractive 
interactions, and recent applications suggest that even more rapid rate decays are also physically 
realistic. 

PACS numbers: 05.40.-a, 05.50.+q, 68.43.Jk, 68.43.De, 81.15.Aa 



I. INTRODUCTION 



Domain growth in far from equilibrium conditions is 
observed in phase separation of mixtures, dynamics of 
glasses and island coarsening duringor after deposition 
of a thin film, among other systems [H, @, H, 0| • This mo- 
tivated the proposal of many statistical models which ex- 
hibit growth laws for the typical domain size in the form 
I ~ t z , where z is a coarsening exponent [3j. For instance, 
when a system is quenched from an homogeneous phase 
into a broken-symmetry phase, two universality classes 
are frequently found, one of them of curvature driven 
(or diffusive) growth [1, [a], with z = 1/2, and the other 
of conserved scalar order parameter 0,111, with z = 1/3. 
However, many model dynamics do not obey detailed bal- 
ance and may lead to domain growth with other power 
law forms or with anomalous coarsening, in which with I 
grows slower than any power of time. A continuous range 
of coarsening exponents may be obtained by tuning a sin- 
gle parameter in models with relatively simple physical 
mechanisms, e. g. single particle exchange between clus- 
ters 0, On the other hand, anomalous coarsening 
is found in certain models that mimic glassy behavior or 
phase separation [n|, [13] (such behaviour is also present 
in models with detailed balance under certain conditions 
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[13j|). A range of coarsening behaviors is also obtained 
experimentally, e. g. in recent works on shaken granular 
systems ((logi) 1 ^ 2 ) fl4| . separation of mixtures of milk 
protein and amylopectin (0.04 < z < 0.2) [lj| and air 
bubbles in foams (0.2 < z < 0.5) Despite the va- 

riety of possible scenarios which were already shown in 
the literature, the study of simple models with normal or 
anomalous coarsening is still important because it may 
reveal the basic microscopic mechanisms that lead to cer- 
tain macroscopic behavior. Such basic studies may also 
help the development of more realistic models for a wide 
range of processes, such as those in Ref. [l7j . 

A class of models in which islands grow via particle 
diffusion, attachment and detachment (Ostwald ripen- 
ing) is very important in surface science because they 
can explain many features of submonolayer or multilayer 
growth [H, [HI [2(3]. Even the one-dimensional models 
are important in this field, both as a first step to un- 
derstand realistic two-dimensional systems and as mod- 
els for growth of elongated islands [2l|, [2^, [Hj]. These 
one-dimensional models may usually be mapped onto 
zero-range processes (ZRP), whose universal and non- 
univeral prop erties were intensively studied in the last 
years [U, (25|. Here, we will analyze the coarsening pro- 
cess in a class of conserved one-dimensional models with 
those mechanisms, as illustrated in Fig. la. The map- 
ping to a column problem, which is a ZRP, is shown in 
Fig. lb. Isolated adatoms diffuse with unit rate and at- 
tachment occurs immediately after a particle reaches the 
border of a cluster. We study here the case in which the 



rate of detachment from a cluster decreases with increas- 
ing cluster size as an inverse power law of the form 
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with k > 0. We consider a very large lattice (infinite 
for practical purposes), where a non-trivial, continuous 
coarsening process is observed if the system begins in a 
completely random configuration. 

This form of detachment rate could apply with some 
type of long-range attraction between the particles in a 
cluster 26] . This mechanism may not be generic for usual 
surface science applications, but the form may neverthe- 
less be a reasonable approximation for a range of cluster 
sizes. Moreover, it may find applications in other fields, 
such as granular systems, where rates with much faster 
decay (7 ~ exp (— m 2 )) were already used to model real 
systems [lij ]. This is an important motivation for this 
study, and additional support to this claim is provided 
by some of its unusual features. First, cluster growth 
shows features that resemble other ZRP with biased dif- 
fusion [27], [HI because there is a preferential flux from 
the small to the large clusters, despite the model rules 
being completely symmetric. The coarsening exponent 
is z = l/(k + 2), but there is a logarithmic correction to 
the dominant power-law coarsening. Thus, as k — ► 0, we 
obtain z — » 1/2, instead of the value 0=1/3 obtained 
with symmetric rules in Ref. [23| (constant 7) and Refs. 
[27l 128} (decreasing 7, but 7(772) — > 1 as to — > 00). On 
the other hand, the logarithmic correction represents the 
crossover from symmetric to effectively assymmetric par- 
ticle flux which occurs during the exchange of particles 
between neighboring clusters. Another interesting fea- 
ture is the difference between the scaling of the average 
cluster size (all clusters) and the scaling of the typical 
size of large clusters for k < 1, due to the presence of 
high densities of small clusters dominating that average. 
This contrasts to related models, including those with 
deposition and/or fragmentation, whose relevant cluster 
sizes are described by a single scaling relation. These fea- 
tures are accompanied by cluster size distribution with 
non-usual features, including a high negative skewness 
near the typical growing size. 

At this point, it is also important to recall the differ- 
ences from previously studied models with similar mech- 
anisms. The case of constant detachment rate (more pre- 
cisely, k = and 70 <C 1) was considered in Refs. [23l.l29|. 
and shows a coarsening with exponent z = 1/3 up to a 
characteristic time of order 7o -5 ^ 2 - Models with 7(771) 
increasing with m were also analyzed in previous work 
[30| and have prospective application to island formation 
in heteroepitaxy, particularly due to the possibility of 
changing the shape of the island size distributions (from 
monotonic to peaked ones) by tuning temperature or cov- 
erage. In those cases, steady states could be attained in 
infinitely large lattices, but the present model (decreas- 
ing 7(m)) shows a steady state only in a finite lattice. 
The properties of this steady state can be exactly pre- 
dicted from a mapping onto a ZRP: for any rate of the 
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FIG. 1: (a) Illustration of the diffusion (m = 1) and de- 
tachment (m > 1) processes of the model, with the asso- 
ciated rates 7(777). (b) Examples of detachment processes 
(1,2) and diffusion processes (3,4) of shaded particles, in the 
original cluster picture and in the corresponding column pic- 
ture. Dashed lines show the correspondence between clus- 
ter+vacancy and a column in the two pictures. 



form in Eq. |T]) , there is condensation into a single cluster 
whose density tends to 1 as the lattice size increases [24| . 

Our results for the average cluster sizes, including the 
logarithmic corrections to the dominant behavior, will 
be derived from a scaling theory presented in Sec. Inland 
will be confirmed by simulation data. In Sec. IIV1 we 
will write the Master equation of the process in an inde- 
pendent interval approximation (II A), and obtain some 
exponent relations. However, because of its neglect of 
important correlations, some results of this IIA do not 
agree with the scaling ones. But, after some adjustment 
it is able to predict the correct dominant coarsening ex- 
ponent. The simulation results for cluster size distribu- 
tions are shown in Sec. |Vl which qualitatively confirm 
the assymmetry predicted by the IIA and the proposed 
scaling relations for small and for typically large clusters. 
Finally, in Sec. IVI1 we present our conclusions. 



II. SCALING THEORY 



A. Basic definitions and coarsening with constant 
detachment rates 



Here we review the heuristic scaling approach based 
on random walk properties used to predict the time evo- 
lution of the typical cluster size. We consider the model 
with small mass-independent detachment rates, i. e. 
7(771) = 70 -C 1 for to > 2, while 7(1) = 1 (free particle 
diffusion). These arguments were formerly presented in 
Ref. [1^1 and follow similar lines of those applied to other 
ZRP in Refs. [HI, HI] . We denote the typical cluster size 
as M, which must be understood as an average over the 
largest (time-increasing) sizes which are statistically rel- 
evant. This average excludes, for instance, clusters with 
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FIG. 2: (a) Configuration of two neighboring large clusters 
with typical size M, separated by an empty region of size I. 
(b) Scheme with a sequence of configurations during the mi- 
gration of a particle from the left to the right cluster, (c) 
Scheme with successive migrations of particles from one clus- 
ter to the other, until the right cluster doubles its mass at 
the expense of the left one. (d) Scheme with the number of 
particles ANl and AN a detached from neighboring clusters 
during a time interval At. From the model rules, the smaller 
cluster (left) loses more particles than the larger one (right). 

Fig. 2a shows two neighboring clusters of size M sep- 
arated by a gap of size I — rM, where r is related to the 
particle density (coverage) 9 by 



(2) 



1 - 



A characteristic time At coa is that in which such clusters 
exchange so many particles that one of them approxi- 
mately doubles its mass at the expense of the other. This 
time is estimated below. 

The time for detachment of a single particle from the 
edge of a cluster is of order r ~ l/ 7 o- However, after 
detachment it is much more probable for this particle to 
reattach to that cluster than to diffuse to the other clus- 
ter. The probability of traveling a distance I before going 
back to the original cluster is l/l, as determined by the 
solution of "the gambler's ruin problem" [3l[ - see also 
Refs. [23|, [2|| . This means that the particle will detach 
and reattach to the original cluster a number of times of 
order / before migrating to the neighboring cluster. This 



is illustrated in Fig. 2b, where for simplicity only two un- 
sucessful detachments (i. e. detachment-reattachment), 
labeled (1) and (2), were shown. Consequently, success- 
ful migration of a single particle from one cluster to the 
other takes place after a time At m i g given by 



Air 



rl ~ rM/70. 



(3) 



The additional time for random walk of the free particle, 
trw ) is negligible during coarsening. 

The above reasoning implies that single particle ex- 
change does not depend on the current size of each clus- 
ter, but only on their separation I, which is kept fixed 
during the process. This symmetric random exchange is 
illustrated in Fig. 2c. After the migration time r m i g , the 
size of each cluster increases or decreases by one unit with 
equal probability. Thus, in order for the size of one of the 
clusters to increase from M to 2M (and the size of the 
other cluster to decrease from M to zero), the exchange 
of nearly Af 2 particles is necessary. Thus, the coarsening 
time is 



At coa ~ M 2 At mig ~ rM 3 /7o- 
This gives a scaling equation 

M 



dM 

dib At coa 



from which we obtain 
M 



1/3 



(4) 



(5) 



(6) 



Notice that the random walk of the free particle be- 
tween the neighboring clusters takes a time of order 



trw -I 2 - (rM)" 



(7) 



If trw •C t, then during the successful migration time 
there will be only one free particle between the clusters, 
as assumed above. Otherwise, if trw ~ t, it is probable 
that two free particles meet, which leads to the formation 
of an intermediate cluster with those particles. Since the 
time necessary for the small intermediate cluster to break 
is of the same order as the detachment rates from the big 
clusters (r ~ I/70), the coarsening process ends. In this 
situation, we have M ~ — K-r* for the average cluster size 



B. Coarsening with decreasing detachment rates 

Here we extend the previous approach to the case of 
decreasing detachment rates (Eq. [I]). In this case, the 
characteristic time for single particle detachment from a 
typical cluster of size M is 



t(M) 



7 (M) 



M k 
7o : 



(8) 
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where we used 70 ~ 1. 

In contrast to the model with constant detachment 
rates (Sec. Ill Ap , here we observe that coarsening will 
not end in an infinitely large lattice because, during the 
exchange of particles between neighboring clusters, the 
time necessary to break the intermediate cluster is of or- 
der 1, which is much smaller than the detachment time. 
In a finite lattice, this leads to condensation of a finite 
fraction of the particles into a single cluster (with the 
present rates, this fraction tends to 1 as the size increases) 

The time for successful migration from one cluster to 
the neighboring one is 



rM fc+ 7 7o 



(9) 



which now depends explicitly on the mass of the cluster 
from which it detached. Detachment from large clus- 
ters is slower, thus there is a preferential flux of parti- 
cles from small to large neighboring clusters. Eq. (j4]) 
is no longer valid because the number of single particle 
exchanges necessary for two clusters to coarsen is much 
smaller than M 2 . When the neighboring clusters have 
nearly the same size, random exchange of particles takes 
place, but as soon as the sizes are unbalanced the net 
flux becomes asymmetric. 

The next step is to calculate the number of exchanged 
particles within a time interval At if the mass is unbal- 
anced by an amount x, as shown in Fig. 2d. The numbers 
of detached particles from the left and the right clusters 
during that time are, respectively, 
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j(M-x) At/ (rM), 



j(M + x)At/(rM). (10) 



Consequently, the mass difference x increases by 



Aa: = AN l— AN 1 



7o 



rM k + x 



\ MJ V MJ 



Mj 

(11) 

within time At. 

The time for a net flux of a fixed mass Ax decreases as 
x increases, which means slow coarsening for clusters of 
nearly the same size and rapid coarsening with one big 
and one small cluster. Transfer of unit mass (Aa; = 1) 
takes place in a time of order 
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and the coarsening time is 



= M 



At coa = / J i±l I 
x=l 



rM k+ 
70 



2 r l 



du 



l/M (1 — u) k - (1 + u) k 

(13) 



For typical masses M ^> 1, the integral in Eq. ((T3|) 

u y k -(l + u)- k « 
' 1, we obtain 



is dominated by u <C 1, where (1 - 
2ku + O (it 3 ) . Since we consider k 



At C oa M fc+2 lnM 

7o 



(14) 



Notice that u <C 1 in Eq. ([13")) , which leads to the loga- 
rithmic correction in Eq. (|14[) . physically corresponds to 
the regime of symmetric particle exchange, i. e. neigh- 
boring clusters with approximately the same size. Simi- 
lar arguments were used to calculate coarsening times in 
Ref. (34]. Since k > 0, we observe that At coa is always 
larger than the time for random walk between the clus- 
ters, given by Eq. (JT]), thus particle detachment is always 
the leading contribution to the coarsening time of large 
clusters. 

Substituting Eq. (TJ3J) in the scaling equation (f5]), we 
obtain 



M 



7o t 
r hit 



1 



(15) 



In order to test these predictions, we performed numer- 
ical simulations of the model for several values of k in the 
range [0.25, 3], with coverages 9 = 0.8, in lattices of sizes 
from L — 8192 to L = 32768, so that finite-size effects are 
negligible. Simulations for some smaller coverages were 
also performed, but the coarsening process usually takes 
place at much longer times. The average cluster size (to) 
was obtained from at least 100 configurations for each K, 
up to times of order t = 10 6 . 

Estimates of the exponent z are usually obtained 
from extrapolation of effective exponents calculated from 
(m) (t) . Without accounting for logarithmic corrections 
in Eq. (fT5|) . we define the effective exponents as 



Z eff,l 



In [(m) (t)/(m) (t-St)} 
In \tj (t - St)} ' 



(16) 



with fixed St. On the other hand, in order to account 
i^r the logarithmic corrections in Eq. (TT5)) . the effective 
exponents must be defined as 



In [(to) (t)/(m) (t-St)] 
\u[{t/\ut) / ((t- St) /In (t-5t))\ 



(17) 



%eff,i is plotted in Fig. 3 function of 1/t for 

fc = 3, k = 2, k = 1 and k = 0.25, and z e //,2 is plotted in 
Fig. 3b for the same values of k. Predicted asymptotic 
values z = l/(fc + 2) (Eq. EJ) are 0.2, 0.25, 0.333 and 
0.444, respectively. For all k > 1, we observe that con- 
vergence to the asymptotic z (as 1/t — ► 0) is faster with 
z eff,2- This justifies the theoretically predicted logarith- 
mic corrections. 

However, for k = 0.25 we observe that both z e ff i and 
z e ff,2 converge to z ss 0.12, which is very far from the 
predicted value of Eq. (fT5|) . In Sec. Ill Cl we will show 
that for k < 1 the coarsening exponent for (to) is actu- 
ally different from z — l/(fc + 2) due to the large density 
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of isolated particles. Thus (m) is very different from M, 
which represents the typical size of large, increasing clus- (a) 
ters. However, we will show that M still coarsens with 
the exponent given by Eq. (jT5j) . 
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FIG. 3: (Color online) Effective exponents z e //,i (a) and 
z eff,2 (b) of the average cluster size (average over all clus- 
ters) as a function of inverse time, with coverage 6 — 0.8: 
k = 3 (squares), k = 2 (up triangles), k = 1 (crosses), and 
k = 0.25 (down triangles). 



C. The role of isolated particles 

The successful detachment of a particle from a clus- 
ter, which allows the migration to the neighboring one, 
takes place after a time interval given by Eq. ©. How- 
ever, this time measures the average residence time of the 
particle attached to the original cluster. The total time 
of migration of a single particle has to include the ran- 
dom walk time between the neighboring clusters, which 
is given by Eq. ([7]). 

If t m ig > trwi then the random walk is rapid, thus 
it is very rare to observe a single free particle between 
any pair of clusters and even rarer to observe two. This 
condition is satisfied when k > 1. Fig. 4a shows some 
snapshots of the simulation for k = 2, which confirm this 
behavior. Thus, the large clusters with mass of order M 
are statistically dominant, i. e. M actually represents 
the average cluster mass among all clusters, which we 
denote by (m). 

On the other hand, if k < 1, a large time is spent in the 
random walk between neighboring clusters. During this 
time, the successful detachment of other particles is pos- 
sible (we recall that intermediate small clusters rapidly 
break for 70 ~ 1). This is illustrated for k = 1/2 in the 
snapshots of Fig. 4b. The number of free particles during 





FIG. 4: (Color online) Sequences of configurations (from top 
to bottom) of a certain region of the lattice for (a) k = 2 
and (b) k — 0.5. In both cases, the coverage is 6 — 0.6 
and snapshots are separated by a time interval 10 (simulation 
times are of order 10 4 ). 



trw in the region between two large clusters is of order 

JVi ~ T RW /t mig ~ M l ~ k (18) 

and the corresponding density of free particles is 

pi - N\/M ~ M~ k . (19) 

However, the density of large clusters, whose typical 
mass is M, varies as 



Plarge ~ 1/M. 



(20) 



This means that the free particles (or small clusters 
formed by their attachment) are statistically dominant 
for k < 1. In this situation, M represents the average 
size of large clusters, but not the average size among all 
clusters, which is (to). 

For k > 1, Eq. (fTS"]) is also valid as a density averaged 
in space and time (during most of the time, there is no 
free particle between the neighboring clusters) , thus large 
clusters of size M are statistically dominant and (to) « 
M. 

These results do not invalidate the arguments of Sec. 
Ill Al for the scaling of M, which is still expected to follow 
Eq. (fT5")) for k < 1. The average cluster size calculated 
among all clusters, including free particles, is obtained 
from an average in the region between two large clusters: 

1 • JVi + M ■ 1 M 

(m) ~- 



Nx + 1 1 + M 1 - 

With k < 1, this global average scales as 

InJ ' ZG ~~k + 2 



(21) 



(fc<l). 
(22) 
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This explains the discrepancies in the numerical esti- 
mates of coarsening exponents for k < 1 (Sec. IIIBI) . For 
instance, for k = 0.25, Eq. (f22f predicts zq = 0.111, 
which is consistent with the trend of the data in Figs. 3a 
and 3b. 

In order to test the predicted scaling of M(t), we cal- 
culated numerically average cluster sizes from contribu- 
tions of large clusters only (masses m > 12 for k = 0.5, 
m > 25 for k — 0.25). Corresponding effective exponents 
are defined as 



In [M (t) /M (t - St)} 



In [(t/ hit) /{{t- St) I In (t -St))]' 



0.5 



(23) 
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FIG. 5: (Color online) Effective exponents z e ff,3 for average 
cluster sizes excluding small clusters, with coverage 9 — 0.8: 
k = 0.5 (asterisks) and k = 0.25 (down triangles). 

Zeff,3 is shown in Fig. 5 as a function of 1/t for k = 
0.25 and k = 0.5. Good agreement with the predicted 
asymptotic value z = 0.4 for k — 0.5 is obtained. For 
k = 0.25, the trend of z e ff,3 as 1 /i — > is not consistent 
with the predicted value 0.444, which is probably due to 
corrections to scaling. In both cases, effective exponents 
not accounting for the logarithmic corrections (similarly 
to z e //,i - Eq- H6p show larger discrepancies from the 
theoretically predicted values of z. 

Additional support to our theoretical predictions is 
provided by the numerical study of the scaling of the 
density of free particles. From Eqs. (JTSJ) and (TT5)) , we 
obtain 



Pi 



t 
hit 



-z(i) 



with 



z(l) 



k + 2 



(24) 



(25) 



(i. e. z (1) = zq for k < 1). In Fig. 6 we show 
[t/hi(t)] z{1) pi versus 1/t for k = 2, k = 1 and k = 1/2, 
using the exponents z(l) given by Eq. (|25|) . The con- 
vergence of that ratio to finite non-zero values as t — > oo 
confirms the expected scaling. 
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FIG. 6: (Color online) Simulation results for the scaled den- 
sity of isolated particles as a function of inverse time, with 
coverage 9 = 0.8: k — 2 (triangles), k = 1 (crosses) and 
k = 0.5 (asterisks). 

The densities of other small clusters can be obtained 
from pi by observing that they have high detachement 
rates and, consequently, they may be viewed as a set 
of nearly free particles at consecutive lattice sites. This 
reasoning gives the density of clusters of size m, for m ~ 1 
as 



with 



z (m) = mz (1) . 



-z(m) 



(26) 



(27) 



Simulations also confirm this result for small clusters, 
such as m — 2 and m — 3, for several values of k. 



III. RELATION TO OTHER MODELS 

Our model may be mapped onto a column problem 
which clearly shows that it is a ZRP. A cluster of length 
m in the original problem and the vacant site at its right 
side is represented by a column of mass m in this new 
picture. The mapping is illustrated in Fig. lb. Sets of n 
consecutive vacancies in the original problem are repre- 
sented by n — 1 vacant columns in the new picture. The 
detachment and diffusion processes correspond to hop- 
ping of a particle from a column to the neighboring one. 
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The mass-dependence of detachment rates is translated 
into mass-dependent hopping rates 7(777) = 2e(m) in or- 
der to account for the detachment in two edges of each 
cluster, each one with rate e(m). 

In a finite lattice, condensation of a finite fraction of 
the mass in a single cluster is expected for all densities if 
7(777) — + for m oo. Moreover, the density of particles 
out of the condensate decreases as L — > oo, as explained 
in Ref. [HI- This is the case of our model, and our 
simulations in small lattices confirm those steady state 
features. 

However, while steady state properties of ZRP can 
be analytically calculated, the coarsening process in in- 
finitely large lattices is much more difficult to predict. 
That is the reason why we use scaling approaches, simu- 
lation and analytical tools based on suitable approxima- 
tions (Sec. HVj) to study coarsening of our model. 

Comparison with related models is interesting at this 
point. Godreche [27[ and Gro/3kinsky et al [28[ ana- 
lyzed the ZRP with hopping rates 7(771) = X + b/m us- 
ing heuristic arguments similar to ours (see also review 
in Ref. [13] )• They considered the cases of symmet- 
ric and asymmetric hopping rates, which lead to average 
cluster size scaling as (to) ~ i 1 / 3 and (to) ~ /j 1 / 2 , re- 
spectively. The symmetric case is somehow equivalent 
to our model with mass-independent detachment rates 
(Sec. Ill Al) . since both have constant and nonzero 7(772) 
for to — ► 00 (very large clusters). 

However, it is important to notice that our model with 
k — > 0, i. e. with very weak mass-dependence of hopping 
rates, has z — > 1/2, in contrast to z = 1/3 which char- 
acterizes constant detachment rates. Both models con- 
sider symmetric hopping rates, but the asymmetric flux 
of mass between the neighboring clusters in our model is 
always present and is responsible for the faster coarsen- 
ing, even if k is very small. In other words, coarsening 
in the model with — > is very different from that with 
k = 0. 

On the other hand, we note that z = 1/3 is obtained 
in our model for k = 1. In this case, the detachment 
rates decreasing with cluster size tend to make the coars- 
ening slower, and balances the effect of the asymmetric 
particle flux between clusters, which favors faster coars- 
ening. For k > 1, mechanisms favoring slow coarsening 
are stronger, thus z < 1/3. For k < 1, mechanisms favor- 
ing fast coarsening are stronger, thus z > 1/3. However, 
both mechanisms are absent in the model with constant 
7 and in the model of Gro/3kinsky et al [28| , both having 
z = 1/3. 

The above discussion leads to the the conclusion that 
the same exponents may be obtained with different mi- 
croscopic dynamics, while apparently similar dynamics 
may lead to very different coarsening exponents. It is 
important that such features are considered if one aims 
to model real systems by ZRP or similar models. 



IV. INDEPENDENT INTERVAL 
APPROXIMATION 

A. General formulation 

The full analytic description of systems with stochas- 
tic processes such as those of our model is provided by 
the Master equation, which is most easily written in the 
column picture of Fig. lb. Previously, this approach was 
used to study the (exact) steady states of related models 
which correspond to ZRP [2^, y(| and the coarsening in 
models with increasing number of particles due to depo- 
sition processes [33j. 

The description of the present model is simplified by 
the fact that the process conserves the total particle num- 
bers N. Thus, using periodic boundary conditions and 
a total number of sites L (lattice length in the original 
cluster picture), and denoting by N t (m) the total num- 
ber of clusters of size to (> 1) at time t, it follows that 
(i) N — J2m=i mNt{m), (ii) the number of spacers in the 
column picture is X)m=i^*( m )' anc ^ ( n ^t( TO ) equals 
the number of columns of size m, for to > 0. Hence, de- 
noting by Nt(0) the number of columns of size zero, we 
have YZ=o N ti m ) = L - N = L{\ - 6) (the last step 
defining the coverage 6 in the original picture). Thus the 
total number of columns (including those of size zero) is a 
constant L c = L — N . The density in the column picture 
is p c = N/L c = r (Eq. |2J. 

The system configuration can be specified by the or- 
dered set of numbers of particles in each of the columns 
in succession: (toi, TO2 . . . to^ c ) = {uii}. The probability 
Ptimi} at time t of the configuration {to^} changes by in 
and out processes. Collecting the effects of all such pro- 
cesses in a time step t — ► t + 1 (see e. g. Refs. [H, l30j l) 
gives the full Master equation 

Pt+i{rrii} - P t {mi} 

1 

$3 [7 ( m i-i + 1) Pt (• • • mi-i + 1, mi - 1 . . .) 
1=1 

+7 (mi+i + l)P t (...mi~ 1, to /+ i + 1 . . .) 

-2 7 [ mi )P t {mi}}6 (mi). (28) 

The theta function above (zero for to < 0, otherwise 
unity) is actually redundant as P t (. . . to — 1 . . .) and 
7(771) vanish for to < 0. 

The Independent Interval Approximation (IIA) as- 
sumes that the configuration probability P t {rrii} can be 

factorised as YlfSi Pt,i (to;). That leads to a reduced 
form of the Master equation in which cluster-cluster cor- 
relations are neglected: 

Pt + i,i (to) - P t .i (to) = At (m + 1, - At (to, , (29) 

where 

A t (m,l) = Pt,i(m) 7 (m)e(m-l) + S m , l7 (l)Pt,i(l) 
- [P t ,i (to - 1) 9 (to - 1) + 5 m>1 P t ,i (0)] J (I) 
= 6{m)At(m,l), (30) 
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with 



A t (m, I) = P t ,i (m) 7 (m) - P t .i (m - 1) J (I) , (31) 
and 

m— 1 

In Eq. (f30|) . the Theta function 9 (to) is zero for to < 0, 
otherwise it is unity. A further reduction results from 
neglecting dependences on the column label I, so P t i (m) 
becomes Pt (to) . This form of IIA gives 

P t+ i (to) - P t (to) = A t (to + 1) - At (to) 6 (to) , (33) 

where 

A t (to) = P t (to) 7 (to) - T t P t (to - 1) , to > 1, (34) 
and 



T t = 22 7 ( m ) -Pt ( m )- 



(35) 



A useful result from the IIA equation (|33f for large 
masses to is 

P t {m')-P t+ i(m')=At(m). (36) 

m'— m 

Hereafter we consider the mass-dependent rates in Eq. 
([I]) for to > 1. Unless otherwise stated, we will proceed 
with developments without dependence on column label 
I, i. e. starting from Eqs. 1(55]). (|3"4"|). (|3"5j) and ([3"6]) . 
with X)m=o Pf ( TO ) = 1 ■ Notice that P t (to) here differs 
from the density p m in Sec. (|TTJ by a constant factor 
1 — 9 due to the different lattice lengths used to normalize 
probabilities in different pictures. 



B. Scaling characteristics 

The late time coarsening of large characteristic masses 
is expected to be described by 



(37) 



The exponents y and z depend on k, and y need not 
equal z because the large masses need not dominate the 
normalisation sums, as shown in Sec. Ill CI The region of 
the cluster size distribution where the scaling equation 
(|3"T)) applies and masses are of order t z is hereafter called 
region S. 

For small to, we have 

P t {m) ~r z{ - m \m<^t\ (38) 

where z(m) is defined consistently with Eq. (I26p . This 
region is hereafter denoted as A. 



Finally, P t (0) may strongly contribute to normalisa- 
tion sums because, as coarsening continues and P(m) at 
small to decreases, Pt(0) will approach 1. So, at late 
times, 



1 - P t (0) - f 



(39) 



which defines za- 



C. Direct results for small clusters 

The IIA equations and the above definitions and prop- 
erties directly lead to some results for small m and large 
times. This is a quasistatic situation in which probabil- 
ities slowly vary in time, thus the left hand side (LHS) 
of Eq. (|3"3"|) is negligible. Since Eq. (f3"3")) is valid for all 
to > this leads to A t (m) ~ 0, and Eq. ([34]) leads to 



P t (m)~P t (0)r t m (TO!) A 



(40) 



Since Pt (0) ~ 1, this yelds the form (pIS)) and confirms the 
relation (|27[) among the coarsening exponents of small to 
given that 



T t oc t 



(41) 



The sizes of the terms on the LHS and on the right 
hand side of Eq. (|33p are respectively, for a given to, 



of order (d/dt) [Pt (to)] 



z(m) _ j-l-mz(l) an j 



A t (m + 1) - t~ z( - m+1 > = £-(™+iMi). The quasistatic 
assumption means that the former is negligible compared 
to the latter quantity at long times, thus 



2(1) <1. 



(42) 



This result is also consistent with the scaling picture of 
Sec. Ulland simulation results. 

The sum in Eq. (|35[) can then be separated into the 
contributions from the two regions, A and S. Eqs. ((3"5)) 
and |gT]) [with ([27))] apply to A and Eq. ((37]) applies to 
S, thus 



TOo(t) 

=1 Jm {i) 
t -z(l) + i _ >+t -b+(*-l)«] 



,-fcj.- 



f / dTO 



(43) 

with 1 < to (t) < t z . Eq. g3]) is consistent with (gT]) if 
z(l) <y+(k-l)z. (44) 



Simulations strongly support Eqs. (|4ip and (|43|) . as 
well as (jlH as an inequality (which is also consistent 
with the scaling theory, as discussed below). This implies 
that the sum in T t is dominated by the small to region 
(actually by just the to = 1 term). The result (I40p . which 



implies P t (TO)/[P t (l)]' 
simulation. 



(to) , is also confirmed by 
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D. Results for large clusters and exponents 
relations 



Here we denote by an d J2s * ne summations with 
respect to m over regions A and S, respectively. 

Consider the sum giving the density in the column pic- 
ture 



E 

m— 1 



mP t (to) 
mP t (to) - 



s 



mP t (to) 



(45) 



(the sums being carried out in same way as those giving 
Eq. I43p . Since p c is constant in time, this is consistent 
with 



y = 2z. 



(46) 



Similarly, the density of clusters in the scaling region 
S is 



s 



(47) 



where we used Eq. (l4"6"]) . 

A further exponent relation follows from Eq. (|39[) and 



1-P t (0) = ^F,(TO)+^P t (TO) = 



.r z = 



which gives 



za — min{z (1) , z}. 



(48) 



(49) 



It turns out that the minumum here is z(l) for k < 1 
and z for fc > 1, where small and large clusters are re- 
spectively dominant (this was shown in Sec. [ill and will 
be confirmed in the context of the II A below). 

These considerations warns us that there are several 
average masses, including 



^ mP t (to)/ ^ p t ( m ) 
s 



t z mp t m/ E Pt ( m ) k <z ' 4 ' 

rn—l rn—1 



(50) 



with za given by Eq. ([49]) . 

Now consider the IIA Master equation in the form Eq. 
(|36j) . and the ansatz for the scaling regime, Eq. ([37]) . 
Replacing the time difference by a derivative and the sum 
over to by an integral, we have (also using Eqs. [Tl 1341 and 
ED 



dm dt 



. TO - fe -...t-^)(i-#- 

OTO 



(51) 



The leading order terms on the RHS cannot cancel, since 
they have different dependences on m, so we can ignore 
the subdominant -J^- (which came from the to — 1 argu- 
ment). With x = m/t z and u = m! /t z , the result is 



du [yf (u) + zuf (u)] = 



. r zk x- k - . . . i~ z(1) t~ y f{x). 



(52) 



The quasi-static results for small to came from achieving 
a cancellation on the RHS. For the large m case, the 
different ^-dependences preclude cancellation, but both 
terms on the RHS have the same dominant order if Eq. 
([21))) is valid. As shown in Sec. [Hi this is consistent with 
our scaling theory and with simulation data. 

However, the dominant ^dependences in the LHS and 
RHS of Eq. ([51 give 



y 



l 



-zk — y 



1 



k+ 1 



(53) 



Comparison with z = ^q^, given by Eq. (|15p and con- 
firmed by simulation, shows that this result is not correct. 
The origin of the discrepancy is an important correlation 
missed by the IIA, as will be discussed below. 



E. Inadequacy of the IIA and a heuristic 
adjustment 

The temporal evolution at large to is being misrepre- 
sented by the IIA because it associates a product weight 
Pt (to) P t (1) to the joint occurrence of a free particle and 
a cluster of mass to, not distinguishing between cases 
where the particle and the cluster are adjacent or well 
separated. These two cases are very different for large 
to because of the small probability of detachment of a 
particle from a large cluster and the high probability of 
the subsequent random walk of the particle finishing with 
absorption at the originating cluster. In the full original 
Master equation [Eq. |28[) . with cluster /column labels 
and without factorisation of probabilities] it is easy to 
identify the random walk steps (through the i-labels, and 
since they occur with rate 7(1))- For comparison, we can 
also see them through A t (to, I) in the IIA version retain- 
ing column labels [Eqs. P9" |) , ([5Tj [l . ([3T [) and ([3"2"[) ]. where 
here the inadequate factorisation has been made (which 
does not properly represent the distortion of the walk by 
the large cluster). 

These problems can be adjusted as follows. The ab- 
sorbing aspect of the random walk of a single particle 
near a large cluster reduces the effective rate of migra- 
tion to another large cluster. Given that their average 
separation increases as their average size, t z , the reduc- 
tion is by an extra factor i~ z , to be introduced into the 
terms on the RHS of Eq. (|5ip (consequently, the RHS 
of Eq. [55] also changes by the extra factor t~ z ). This is 
equivalent to the effect included in the scaling arguments 
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of Sec. HU The consequence is that in place of Eq. (f53|) . 
the power counting gives 



z — y — 1 = — zk — y — z 



k + 2 



(54) 



The extra factors t~ z do not modify the quasi-static form 
for the distribution function Pt(m) at small to, thus its 
introduction is still consistent with Eq. (|23|) . 

Thus, using these heuristic arguments we are able to 
predict the correct coarsening exponent and preserve sev- 
eral exponents relations. However, the changes are still 
unable to predict the logarithmic corrections shown in 
Sec. Ill BL which are related to a crossover from symmet- 
ric to asymmetric particle exchange between neighboring 
clusters. 



F. Cluster size distributions 

Using Eqs. ([25]) and (gT|), the quasistatic result ((4T)|) 
for small masses can be rewritten as 



(55) 



Comparing with t~~ v / (p-) (Eq. [37]) . it can be estimated 
that the crossover between the forms for the region A 
and the scaling region S occurs at to = mo(t) where 



m (t) ~ ~t z [l + 0(t- z \nt) 



(56) 



The form ([55]) first decreases with to (due to the in- 
creasing power of t~ z ) but then turns over into an in- 
creasing function when to exceeds 0(t z ). The mini- 
mum is at to — m(t) such that = [mP(m)] ~ 
[mk (Inm — zhxt — 1)] = k (In to — xlni), so 



m(t) = t z 



(57) 



Thus the minimum is near the crossover region. Simu- 
lations consistently show that the scaling starts just be- 
yond the minimum and that the quasistatic results (I40|) 
and (|55[) work well up to just beyond the minimum. 

Eqs. (|15[) and (|57[) imply that the position of the min- 
imum decreaases with increasing k. This is also seen in 
simulations, and is consistent with small clusters having 
largely to = 1 for k > 1 and a greater spread for k < 1. 

The adjusted form of Eq. (|5"2"|) for the scaling function 
(Sec. ITVEl is, using Eqs. (25) and (06]), 



du [2/ (u) + uf («)] = (ax~ k -b)f (x) (58) 



where the factors of t have consistently cancelled by us- 
ing the correct coarsening exponent (Eq. I54p . and a 
and b are constants associated with 7 (to) and T, respec- 
tively. Differentiating Eq. (f58"| with respect to x gives 

din /Or) = 2-akx~^ +1 '> 
dx x+ax~ k —b ' 



hence 



/ (x) oc exp 



dx 



2 - akx-( k+ V 



(59) 



For small x, the integrand in the indefinite integral 
is dominated by —kx~( k+1 ' ) /x~ k , which integrates to 
lnx~ k , thus 

/(x)oci fc . (60) 

For large x, the dominant part of the integrand is 2/x, 
giving 

/(x)ocx^ 2 . (61) 



V. SIMULATION RESULTS FOR CLUSTER 
SIZE DISTRIBUTIONS 

Despite the problems of the IIA to predict the coarsen- 
ing exponents and the absence of the logarithmic correc- 
tions in the time scaling, even after suitable adjustment 
(Sec. HVEp , it progresses beyond the previous scaling 
theory (Sec. |TT| by providing information on the cluster 
size distributions, which can now be compared to simu- 
lation data. 

The unusual shape of the cluster size distribution in 
this problem is illustrated in Fig. 7 for k — 1 (t = 5 x 10 5 ) 
and k = 0.5 (t — 10 6 ). There is a rapid (power-law) 
decrease of P(m) for small to, usually until to of order 
10, and a peak appears at large to, i. e. in the range of 
typical large clusters. For k > 1, the statistical weight of 
the small clusters decreases with time, i. e. the left side of 
the curve becomes smaller when compared to the peaked 
region. For k < 1 the opposite occurs: as time increases, 
the weight of the small to region increases and the peak 
becomes relatively smaller. Indeed, the curve for k — 1/2 
in Fig. 7 shows that the probability of isolated particles 
or dimers is 100 to 1000 times larger than the probability 
of sizes in the peaked region [for instance, P(l) « 0.77]. 

The first important result of the IIA is Eq. (|57|) for 
the scaling region (the region of the peak in Fig. 7), 
with y given by Eq. (|46[) . Simulations show that this 
result is valid with t replaced by t\ = tj lni, which is an 
expected correction. This is illustrated in Figs. 8a and 
8b, where we show log [tfPt (to)] as a function of m/ti z 
for k — 0.5 and k — 2, respectively, and three different 
times for each k. The good data collapse (particularly 
for the largest times) is obtained with z = 1/ (k + 2) and 
y = 2z, as predicted by Eqs. (fT5)) and (|4*o| . 

A power-law in the left tail of the scaling function f(x) 
is observed in our simulations, but the exponents are dif- 
ferent from those predicted in Eq. (|60|) for small k. For 
instance, for k — 0.5, the exponent is 1.07. For larger k, 
the agreement is slightly better, e. g. exponent 2.05 for 
k = 2. Anyway, one interesting feature of the IIA results 
(|60"f and (|6"Tj) is that the left tails of the distributions are 
heavier than their right tails for k < 2. In other words, 
the distributions have negative skewness. This is clearly 
observed in Fig. 8a, for k = 0.5, while for k = 2 (Fig. 
8b) the skewness is closer to zero (but still negative). 

The negative skewness of cluster size distributions is an 
uncommon feature in this type of problem in one dimen- 
sion; for instance, the distributions in coarsening with 
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FIG. 7: (Color online) Cluster size distributions for k 
at t = 5 x 10 e 
curve) . 



(dashed curve) and k = 0.5 at t = 10 (solid 



O 




FIG. 8: (Color online) Scaled cluster size distributions in the 
scaling region for: (a) k = 0.5 at t = 10 (crosses), t — 2 x 10 
(squares) and t — 10 6 (triangles); (b) k — 2 at t = 2 x 10 5 
(squares), t = 10 6 (crosses) and t — 5 x 10 6 (triangles). 



constant detachment rates are positively skewed [29j |. as 
well as those in the steady states with some rate func- 
tions which increase with cluster size (due e. g. to re- 
pulsive interactions) [3(3]. Thus, in a real system, that 
feature would suggest the presence of attractive interac- 
tions leading to a decrease of the detachement rate with 
cluster size. On the other hand, it is important to notice 
that it is a common feature in two dimensions, both in 
point islands models (which are two-dimensional ZRP) 
and in extended islands models [20| . 

The scaling of small masses (Eq. I55[) is confirmed in 



Figs. 9a and 9b for the same values of k, again with 
the logarithmic corrections in the time t. There we plot 
log [Pt {rn)\ versus m log (m/t z ), which is proportional to 
the logarithm of the RHS of Eq. (|55|) . In Figs. 9a and 
9b, one important point is the large range of both vari- 
ables (horizontal and vertical), which span 2 to 6 orders 
of magnitude. This explains the discrepancies from a per- 
fect data collapse when compared to the scaling regime 
in Figs. 8a and 8b. 



6, 

O 
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(a) 



-3 -2 -1 
mlog(m/y) 
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FIG. 9: (Color online) Scaled cluster size distributions in the 
small mass region for(a) k = 0.5 and (b) k = 2. Symbols are 
the same of Figs. 9a and 9b. 



VI. CONCLUSION 

We studied conserved one-dimensional models of par- 
ticle diffusion, attachment and detachment from clusters, 
where the detachment rates decrease with increasing clus- 
ter size as j(m) ~ m~ k . Heuristic scaling arguments 
based on random walk properties were used to predict 
the scaling of the typical cluster size as (t/\nt) z , with 
z = 1/ (k + 2). The coarsening of neighboring clusters 
is characterized by initial symmetric flux of particles be- 
tween them followed by an effectively assymmetric flux 
due to the unbalanced detachment rates (despite the 
symmetric model rules). For k < 1, the average cluster 
size does not scale as the size of typically large clusters 
due to the high densities of small clusters, which domi- 
nate that average. We also solve the Master equation of 
the model under an independent interval approximation, 
which predicts some exponent relations and the correct 
dominant coarsening exponent after suitable changes to 
incorporate effects of correlations. These results are con- 
firmed by simulation, which also shows the negatively 
skewed cluster size distributions (particularly for small 
k) and the different scaling relations followed by small 
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clusters (sizes of order 1) and by typically large clusters 
(size of order t z ). 

The rate functions analyzed here may arise from as- 
sociating (Arrhenius) detachment rates with potentials 
U(m) for particles at the end of a cluster of size m. U (m) 
is then a sum, from I = 1 to m — 1, of pair potentials V(l) 
for separation I with Coulomb-like (inverse of distance) 
attractive form. In a real system, such interaction is not 
expected to be valid for all sizes, but may be a reason- 
able approximation for some ranges, in a similar way that 
long range repulsion between adatoms on a surface rep- 
resents substrate-mediated interactions. The particular 
coarsening features discussed here will certainly help to 
identify such application. 

It is also interesting to note that our model with 
7(771) ~ exp (— m 2 ) was already studied in Ref. [HI and 
quantitatively describes experiments with a shaken " gas" 
of steel beads distributed among a set of boxes. The av- 

1 /2 

erage cluster size increases as (logt) and the density 
of particles in the boxes without big clusters decrease as 
1/t. These results can be obtained by a direct extension 
of the scaling arguments of Sec. [II] (the second one may 
be viewed as the k — > 00 limit of Eqs. [Ml and l25l) . 

From the theoretical point of view, this work contains 
some important advances. First, we show how random 
walk properties and simple model rules are able to pre- 
dict the coarsening law including a logarithmic correc- 



tion, which is a non-trivial task at the level of a scaling 
theory. Moreover, this correction is shown to be a conse- 
quence of a continuous competition between symmetric 
particle flux between neighboring clusters and a domi- 
nant assymmetric flux, despite the absence of a spatial 
bias in the model rules, in contrast with other ZRP where 
asymmetric flux appeared only as a consequence of such 
bias. Finally, the different scaling relations obeyed by 
small clusters and by typically large clusters, which en- 
able the former to be statistically dominant when k < 1, 
contrasts with other models with similar physical mech- 
anisms (even those involving deposition and/or fragmen- 
tation), where a single scaling relation is sufficient to rep- 
resent all relevant cluster sizes. 
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